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ABSTRACT 



1 INTRODUCTION 



According to the concordance ACDM model (|Carroll et al. I 
1 19921 ). dark energy (in the form of a cosmological con- 
stant or a varying scalar field) and dark matter (hereafter, 
DM) are the dominant actors on the co s mological scene 
dKomatsu et al. II2OO9I : iPercival et al. II2009I : iLampteil et al. I 
I2OO9I ). In particular DM represents most of the total mass on 
galactic and cluster scales and drives the formation and evo- 
lution of cosmic structures. Roughly speaking, DM haloes 
form when the expanding matter within (and surrounding) 
an overdense region experiences deceleration because of the 
gravitational force, decouples from the Hubble flow, col- 
lapses, and eventually, virializes. Models including all these 
processes have been realized using different approaches. 

N - body simulations have been the primary instrument 
to fully implement the nonlinearities of the formation pro- 
cess which are realized in the dark halo growth. Despite they 
do not allow to catch the full physics of the galaxy formation, 
coUisionless simulations have been successful in reproducing 
a wide range of galaxy properties, e.g. the spherically av- 
eraged halo density profile, pDM(r), which has been found 
to be well described by a double power -law relation with 



Scaling relations among dark matter (DM) and stellar quantities are a valuable 
tool to constrain formation scenarios and the evolution of galactic structures. However, 
most of the DM properties are actually not directly measured, but derived through 
model dependent mass mapping procedures. It is therefore crucial to adopt theoreti- 
cally and observationally well founded models. We use here an updated version of the 
secondary infall model (SIM) to predict the halo density profile, taking into account 
the effects of angular momentum, dissipative friction and baryons collapse. The result- 
ing family of halo profiles depends on one parameter only, the virial mass, and nicely 
fits the projected mass and aperture velocity dispersion of a sample of intermediate 
redshift lens galaxies. We derive DM related quantities (namely the column density 
and the Newtonian acceleration) and investigate their correlations with stellar mass, 
luminosity, effective radius and virial mass. 

Key words: dark matter - galaxies : kinematic and dynamics - galaxies : elliptical 
and lenticulars, CD - galaxies : formation 



Pdm oc r~ in the outer regions and poM oc r~" at their 
centers with the exact value of a rema ining a matter of con - 
troversy. In the popular NFW model (|Navarro et al. |[l997l ) 
they find a — 1 i ndependently on ha l o mass, while either 
a steeper a = 1.5 jMoore et al. 11X99^ : iGhigna et 811120001 : 
Fukushige fc Making I I2OQII) or even shallower values (e.g.. 



Power et al. 1 120031 : iFukushige et al~l |2004 iNavarro et al. 1 
2004 ) have been claimed elsewhere. It is also possible that 
Q is not universal at all, but rather depending on halo 
mass, merger histor y and substructures (|jing fc Suto |[2000l : 
iKlvpin et aTlbOOll ). 

On the contrary, semi - analytical models are more fiex- 
ible offering the possibility to i nclude a vast va r iety o f phys- 
ical in gredi ents. In partic ular, iGunn fc Gott I (|l972l ). [Gottl 
(|l975l ) and iGunnI (| 19771 ) introduced the secondary infall 
model (SIM) to describe the collapse and virialization of 
halos that are spherically symmetric, have suffered no ma- 
jor mergers, and have undergone quiescent accretion. After 
these first analysis, other works have relaxed the assumption 
of purely radial self - similar collapse by including non-radial 
motions arising from secondary perturbations and taking 
care of both angular momentum and stars to lead shallower 



2 V.F. Cardone et al. 



or st eeper density profile s depending on the halo mass (see, 
e.g.. iDel Popolo (|2010ah and refs. therein). 

Numerical and semi - analytic models generally lead to 
different predictions which need a detailed observational 
scrutiny. Furthermore, the adoption of the proper density 
model is the basic ingredient to derive the global DM prop- 
erties which are critical parameters in the galaxy formation 
scenario; e.g., the virial mass which is considered the driver 
of the h eating process that might affect the star form ation 
history (jPekel fc BirnboimI 12009 : ICattaneo eraDl2008l') and 
the overall star formation efficiency ( Conrov &l WechsleJ 
I2OO9D . 

From this point of view, scaling relations among DM 
and stellar quantities may provide an important test to con- 
strain both the formation scenarios and the DM properties. 
Early -type galaxies (ETGs) are ideal tools for these aims. 
First of all they are found to lie on the so-called funda- 
mental plane (FP) which tightly relates the central velocity 
dispersion, effective radius and surface brightness. In par- 
ticular, the well known deviation (or "tilt") of the FP with 
respect to the expectation of the virial theorem is still to 
be fully understood. The different explanations proposed 
rely in turn on non homology, variation of the stellar M/ L 
ratio with luminosity and varying DM content (see, e.g., 
Busarello et al. 1119971. iD'Onofrio et al.ll2006l , iTortora et 



20091 . T-l-09, hereafter). Each of these solutions may tell a 



different story about the interplay between the DM and the 
stellar component, thus, it is clear how the constraints on 
scaling relations can help to shed light on formation and 
evolutionary processes. 

Recently, there have been growing evidences that 
DM mass fraction in the ETGs central regions is 
an increasing function of stellar mass (or luminosity) 
hence supporting the idea th at DM mig ht be t he main 
drive r of the FP tih ICappcllari et all [2006._ Bo ltoiTet al.l 

■ M 



20071. iHvde fc Bernardir[2009l. ICardone et al. I (20091. T-H09 



Auger et al.ll2010l '). On the other hand, the mean 3D DM 
central density has been found to decre ase with mass 
and luminosity (|Thomas et al.ll2009l . T+09. iTortora et al. I 
I2OIOI ). while there are contradictory results on the uni- 
versality of the column density Sum = MDM,proj/7rii^ 
(where MoM.proj is the projected DM mass within the ra- 
dius R) with some results arguing fo r its constancy over 
12 orders of magnit ude in luminosity (|Donato et al. I [20091 : 
iGentile et £11120091) and other works fi nding a correlation 
with halo mass iBovarskv et al. I (|2009l '). Part of this con- 
troversy may probably be ascribed to the different assump- 
tions on the halo model and stellar initial mass function 
(IMF) or t he adopted scale radius, a s recently argued by 
some of us (jCardone fc Tortora II2OI0I . hereafter CTIO). On 
the other hand, it is also possible that Sum changes with 
the morphologica l type as s uggested by the recent results 
in [Napolitano. R omanowskv fc Tortora (2010j), where the 
central projected density in ETGs is found to be, on aver- 
age, systematically higher than the same quantity for spiral 
and dwarf galaxies (see also B09). 

In order to further investigate this issue, we present 
here the analysis of the above scaling relations based on the 
SIM density profile obtained in Del Popolo (2010) adding to 
the usual recipe of the gravitational collapse, the effects of 
ordered and random angular momentum, dynamical friction 
and adiabatic contraction due to the baryonic collapse. 



We use Einstein radius and velocity dispersion data 
from a sample of intermedi ate redshift (Iz) .2) lens galax- 
ies from SLAGS survey (|Auger et al. 1 12009| . A-l-09, here- 
after) to constrain the model parameters and derive different 
scaling relations. A general overview of the model is given 
in Sect. 2, while in Sect. 3 we introduce the lens sample and 
describe the fitting procedure. Our main results are shown 
in Sect. 4 and discussed in the concluding Sect. 5. 



2 THE HALO MODEL 

The density profile of DM haloes is here obtained by us- 
ing the analytical method described in Del Popolo (2009, 
hereafter DP09) which we refer the interested reader to for 
more details. Here, we give a brief descriptions of the model 
properties which are of main interest to our aims. 

The hal o profiles are derived by assuming the secondary 
infall model (|Gunn fc Gott iri972l ) where a bound mass shell 
of initial comoving radius Xi expands up to a maximum ra- 
dius (or turnaround radius) Xta. As successive shells expand, 
they acquire angular momentum and then contract on orbits 
determined by the angular momentum itself, while dissipa- 
tive processes and eventual violent relaxation intervene to 
virialize the system converting kinetic energy into random 
motions. The final density profile may then be computed as : 



p{x) = 



(a;/a;ta)^ 



1 + 



d\n{x/xta 
dlnxta 



(1) 



with /Ota(2:ta) the density at turnaround and x/xt^ referred to 
as the collapse factor (see Eq. A18 in DP09). To describe the 
proto-haloes density profile, DP09 considered the profile of a 
peak in the density field generated according to the Baarden 
et al. (1986) power spectrum and then took into account 
angular momentum, dynamical friction and the presence of 
baryons following the steps described below. 

First, the angular momentum is decomposed in an or- 
dered component, related to the tidal torques experienced 
by proto-haloe s, and a random compo nent connected to ran- 
dom velocities (|Rvden fc Gunnll 19871 ). The ordered term is 
computed following Ryden (1988), while the random part is 
assigned to proto-structures according to Avila- Reese et al. 
(1998). A term related to the dynamical friction force has 
been expl icitly introduced in the equations of motion and 
evaluated iKandrup I (|198G| ) dividing the gravitational force 
into an average and a random component generated by the 
clumps in the hierarchical universe. Finally, some adiabatic 
contraction of the halo, due to the baryonic collapse, has 
been taken into account through the formalism of Klypin 
et al. (2002) and Gnedin et al. (2004), also including the 
exchange of angular momentum among baryons and dark 
matter. 

The final product of this halo formation method gives 
the DM density profile as function of the radius r and the to- 
tal halo mass Mvir. The latter is the only parameter needed 
in order to specify the halo density, being the halo inner 
slope a a function of the virial mass as well. As discussed 
in DPIO, the dependence of a on Mvir breaks the univer- 
sality of the halo profiles and favour Burkert (1995) models 
at dwarf scales, and models steeper than NEW at normal 
galaxy scales. 

For the analysis we want to propose in the following. 
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it is more convenient to handle some analytical halo den- 
sity profile, thus we decided to approximate the numeri- 
cal DPIO models w ith a generalized NFW density profile 
l|jing fc Suto llioocj) which allows to accommodate varying 
inner slope a beingj: 

where pcrit{z) = 3H^ (z)/8tvG is the critical densit^Q of the 
Universe at redshift z and Cvir ~ Rvii/Rs is the halo concen- 
tration (with Rs and -Rvir the radius where the logarithmic 
density slope equals -2 and the virial radius respectively). 
Differently from the original SIM we want to approximate 
(which is fully assigned by the halo mass only), the general- 
ized NFW profile is formally a function of two parameters, 
namely (cvir, A/vir)- Thus, in Eq.((2]), the dependence on Cvir 
is fictitious and we can fit the same halo profile with dif- 
ferent Cvir values by changing the corresponding Avir- Due 
to this liberty in the Cvir choice, we therefore arbitrarily 
scale Cvir with the tot al mass Afvir using the popular relation 
l|Bullock et al. II2OOII ) : 



= 12.81 



Mvir 



1012 Mq 



and then fit for Avir as function of Cvir obtaining : 



6.626 X 10" 

Ton 



|c^. 



In (1 -I- Cvir 



1 



(3) 



•(4) 



With the chosen setup, the analytic expression of the in- 
ner logarithmic slope 7 as a function of the mass (Afi2 = 
Afvir/10"'^'^ Mq) turns out to be written as: 



7 •■ 



1.166M, 



1/3 



1 



(5) 



We have checked that this analytical model fits extremely 
well the numerical density profile over the mass range 10^" ^ 
M/Mq ^ 10^* and is fully described by the virial mass as 
a single parameter, as prescribed by the SIM numerical re- 
sults. Moreover, it is worth stressing that the use of the 
above Cvir - M^iv relation is just a convenient choice to sim- 
plify the search for an analytical approximation which does 
not affect the final accuracy of the numerical SIM profile 
fitting 

As a final remark, we warn the reader that the above 
fitting formulae for Avir vs Cvir and 7 vs A/vir have been 
obtained by considering haloes at 2: = (since there is a 
larger statistics and better resolution), while we will adopt 
it also at the intermediate redshifts of the lenses we will 
consider later. As can be seen from Fig. 2 in Del Popolo 



1 The adoption of an analytical approximation to the numerical 
output of the model will make the adoption of the DP09 results 
easier to handle in the mass mapping proposed below. 

2 We assume a concordance ACDM cosmological model so that 
H'^{z)/Hl = 0^(1 + 2)^ + (1 - ^m) with (Qm, h) = (0.3, 0.7). 
^ In principle, one could have adopted a whatever functional form 
provided the relation Avir - Cvir is adjusted in such a way that the 
approximated density profile still fits the numerical one. For this 
reason, one has not to update the Cvir - A/vir relation to account 
for a different cosmology or redshift. In particular, should one 
consider systems with z > 0, one must still use Eq.||3ll without 
scaling Cvir by (1 -|- z)~^ as usually done in literature. 



(2010), the evolution of 7 with z is actually quite small from 
z — Oto z = 1 over the mass range of interest here so that we 
prefer to rely on these well checked approximations rather 
than trying to fit less numerically accurate higher z profiles. 



3 TESTING THE SIM MODEL 

There are two main characteristics of SIM halo model ob- 
tained above that make it particularly interesting : i.) it is 
theoretically well founded and intuitively incorporates most 
of the dark and baryonic collapse physics, and ii.) it ends 
up with a halo family which depends on a single parameter 
(the virial mass) and can be written analytically as a gen- 
eralized NFW (once specified the dependence of the halo 
normalization on the NFW parameter Cvir). 

As a first observational test for this halo model, we 
start with a sample of ETGs for which we can use a multi- 
technique approach as in CTIO. The sample includes 59 
ETGs from the lenses catalog collected by the Sloan Lens 
ACS (SLAGS) survey (Auger et al. 2009, hereafter A+09) 
for which the velocity dispersion dap (within a circular aper- 
ture of radius i?ap = 1.5") and the Einstein radius Re, and 
hence the projected mass within it, Me = Afproj(i?£;), are 
measured. Following A+09, we will model the light distribu- 
tion with a de Vaucouleurs (1948) profile with the effective 
radius Rb and total luminosity Ly set to the values inferred 
from the V - band photometry. Finally, we use the estimate 
of the total stella r mass, Af ^. , from the SLAGS team (A-f 09, 
Table 4) where a lSalpeterl (|l955l ') IMF is assumed. 

The median values of R^p/ Rc and Re/ Rc (respectively 
0.62 and 0.51) indicate that the data probe the galaxy inner 
regions where we need to carefully account for the stellar 
contribution to the model esti mates of (Tap and Me- To this 
end, we adopt the PS model jPrugniel fc Simien Ill997l ) as 
its projection closely mimics the Sersic (1968) surface bright- 
ness profile and provide an analytical form of the physical 
quantities of interest. 

Being the projected profile a Sersic model, the lensing 
properties of the PS model are also analyti cally computed 
iGardone ll20o3 : lEhasdottir fc MoUer llioOTi ). In particular, 
the projected mass within ^ = R/Re is : 



r(b„c 



r(2n) 



(6) 



with Mi, the total stellar mass, T{x,y) and T{x) the incom- 
plete and complete F functions and we set n = 4 to mimic 
the deprojected de Vaucouleurs (1948) profile used by the 
SLAGS team to fit the surface brightness of their galaxies. 

The stellar and DM mass models are finally used as in- 
put for the computation of the luminosity weighted velocity 
dispersion profile. As a first step, the line of sight velocity 
dispersion is given by (|Mamon fc Lokas II2005I ) : 



I{R)aLiR) 



K{v/0p4v)Mtot{v) 
V 



dr] 



(7) 



where I{R) is the Sersic intensity profile, 77 = r/Rc, Mc 
and are the total mass and the stellar density at 
the stellar AI/L ratio, Aftot(??) is the total mass, K{ri/(^) a 
kernel function depending on the choice of the anisotropy 
profile and the tilted quantities are normalized with respect 
to their values at Rc. We consider only isotropic models and 
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Table 1. Main stellar and DM quantities of the lens ETG sample adopted. We first show lenses of the B sample, then the ones to be 
added to get the G and F samples. Stellar parameters are from the SLAGS collaboration (A-|-09), while DM quantities are derived by 
our best fit SIM model to the data. Columns are as follows: 1. galaxy ID, 2. V-band luminosity Ly, 3. logarithm of Re, 4. logarithm of 
the total stellar mass as taken from A-l-09 (assuming a Salpeter IMF), 5. model virial mass, 6. 3D DM mass fraction, 7. logarithm of 
the DM column density iSdm(^g)i 8. logarithm of DM acceleration, 9. logarithm of the Newtonian stellar acceleration. All quantities are 
given with their Icr errors except log Ly and log Rc that have negligible uncertainties. 



take the corresponding K{r]/S,) from Appendix B of Mamon 
& Lokas (2005). The observed quantity is then obtained by 
luminosity weighting aios in a circular aperture of radius 
Ra.p. Note that, according to the SDSS survey strategy, i?ap 
is fixed to 1.5 arcsec so that ^ap = Rs.p/Rc changes from one 
lens to another. 



Having set the stellar component quantities from pho- 
tometry and mass estimates, we are left with only one un- 
known parameter, namely the halo mass Mvir- We find its 
fiducial value minimizing the merit function : 
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where, as detailed in Cardone et al. (2009), the uncertain- 
ties (ffo-, ee) are obtained by summing in quadrature the ob- 
servational errors and the theoretical ones as derived from 
the propagation of uncertainties on (M*, Ra). For each lens, 
the best fit value is the one minimizing X^(-^vir), while 1 
(2) a confidence levels are obtained solving Ax'^(Mvir) ~ 

x'(Mvi.) - Xmin = 1.0 (4.0). 

In order to select only the lenses with the higher signif- 
icance of the model fit, one can roughly set a threshold on 
the reduced values. For our best-fit models, these are typ- 
ically smaller than 1 (for a single degree of freedom - two 
observational constraints vs one model parameter) mainly 
because the inclusion of theoretical errors on (£^,££0. We 
have therefore decided to be conservative and define to split 
the sample on the basis of their observed uncertainties rather 
than the significance of the fit and define a full (F), good (G), 
best (B) as follow. 

The sample F contains all the 59 lenses. The G sample 
is made out of 51 lenses with halo mass in the range 11 ^ 
logMvir 14 and best fit values for ((Tap, Mb) satisfying 



-2.0 SC (a: 



< 2.0 



-2.0 s; {M^''" - A4-f)/£E s; 2.0 . 

The B sample is the most conservative one being made out of 
46 lenses with mass in the same range, but following quality 
parameters 

-1.0 ^{a:^'^ all) 1.0 



-1.0 SC {Me" 



a/e-^)/£e s; 1.0 



The overall agreement of the best fit values of ((Tap, Me) 
with the observed ones is always quite good with rms{l — 
~ 7% (5%) and rms{l- M^^ /Mg'") ~ 13% (9%) 
for sample F (B). As a consequence, most of the results we 
will discuss in the following are quantitatively (within the 
errors) independent on the sample adopted, thus we will 
refer to sample G as the trade-off between fit quality and 
improved statistics. 

Looking more closely to the galaxy filling the F sample 
(i.e. the bottom rows of Table [T]), they mainly differ from 
the G sample for their higher logMvir, exceeding the im- 
posed limit of 1O^*M0 fulfilled by samples B and G. This 
has been set as an upper limit for virial masses of typical 
galaxy systems, and also to match the range of validity of 
our analytical approximation of the SIM (see Sect. 2). 

Such a large log Mvir model estimates might be a warn- 
ing about the accuracy of the approximation adopted to 
convert the SIM density profile in the one parameter gen- 
eralized NFW which tend to overestimate the virial masses. 



* It is worth noting that such a noise term can be hardly reduced. 
To understand why, let us consider the case of the projected mass 
which is the sum of the stellar and DM terms. The stellar term 
is proportional to the total stellar mass and hence is known with 
an uncertainty obtained propagating those on the luminosity and 
the stellar M/L ratio. This later is actually the major source of 
uncertainty in the final budget and cannot be reduced unless one 
re-derive the estimate of relying on a larger set of colours than 
the one used by the SLACS collaboration. 



We have checked though that all (but one) of the seven most 
massive systems filling the F sample reside in overdense re- 
gions (following the definition adopted in Treu et al. 2009) 
so that we cannot exclude that this mass excess in the final 
virial mass estimate is real and given by the contribution 
of the cluster itself which we can not constrain with our 
limited data and is beyond the purpose of this analysis. Fi- 
nally, these systems turns out to have unreasonably large 
virial masses even if we would use the NFW rather than the 
SIM model (see later). Thus, we can confidently conclude 
that they might not represent a critical issue specifically for 
the SIM model, but rather that the lensing model is missing 
some external field to be tracked with forthcoming more ac- 
curate analyses. To be conservative, in the following, we will 
exclude these very massive systems in all the scaling relation 
estimates. 

Since we have only two observed quantities for each lens, 
we expect to determine Mvir on a case -by -case basis with 
some large uncertainty (see e.g. Table [T]). In fact, accord- 
ing to the standard propagation of errors, most of the mass 
related quantities (such as the DM mass fraction and the 
column densities) will be known with poor precision. As a 
possible way out, one could try to improve the constraints by 
binning the galaxies according to luminosity or stellar mass 
and then fitting for a parameter which is as sumed to be the 
same for all the objects in the same bin (e.g.. lCardone et al. I 
I2OO9I ). Unfortunately, this is not possible here because of the 
non universality of the DM mass profile. Indeed, since Mvir 
is obviously different from one lens to another, the slope 7 
and the concentration Cvir will differ too so that the objects 
in the same bin would have intrinsically different proper- 
ties. As a consequence, we do not attempt any binning and 
prefer to deal with large error bars rather than introducing 
possible systematics in the analysis. 



3.1 SIM vs NFW and Burkert models 

The above analysis has demonstrated that the proposed SIM 
model is able to fit the observed aperture velocity disper- 
sion and projected mass within the Einstein radius for the 
SLACS lenses sample with a reasonable significance. One 
might want to check whether other "standard" halo den- 
sity recipes can do a comparable job. We have commented 
in Sect. 2 that our SIM model, because of the dependence 
of the central slope on the virial mass, is able to match a 
wide range of core behaviour from "cuspy" NFW to "cored" 
Burkert models. It is then interesting to compare the best-fit 
results to ((Tap, Me) as performed in Sect. [3] and see whether 
data are able to favor one model with respect the other. 

We need to first remark that the main difference among 
these three models is that NFW and Burkert density profiles 
are intrinsically bi-parametric (although there is a correla- 
tion between the two parameters in both cases, see Sect. 
I3.1.2|) while SIM is mono-parametric by definition. Accord- 
ing with the Ockham's razor principle, this is an important 
argument to take into account in our final considerations. 



3.1.1 SIM vs two-parameter NFW and Burkert profiles 

The two parameters of the NFW and Burkert profiles can 
be similarly recast in terms of the virial mass Afvir and the 
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concentration [fl c = Rvir/Ra where Rg is the radius where 
the slope of the logarithmic density equals -2 in the NFW 
case and the core radius in the Burkert one. Since we have 
only two observed quantities for each lens, it is not possible 
to determine both (c, Mvir) on a case - by - case basis so that 
we must rely on a different strategy. We use the results from 
CTIO, where we have therefore binned the SLAGS lenses 
in ten luminosity bins and expressed (c, Mvir) in terms of 
quantities that can be assumed to be equal for all lenses 
in the same bin. The parameters (c, Mvir) are then deter- 
mined a posteriori for each single lens on the basis of the 
minimization adopted in Sect.O 

Both the NFW and Burkert models fit well the lens 
data, although, as discussed in CTIO, the Burkert model 
leads to quite small virial masses and hence unexpectedly 
low virial M/L values. In order to make a meaningful com- 
parison, we use the same selection criteria defined above 
to define good (G) and best (B) samples for the NFW and 
Burkert fits. It turns out that the number of lenses in the 
G and B samples is (52, 50, 56) and (46, 25, 35) for the SIM, 
NFW and Burkert models, respectively. Note that the same 
lenses enter the G samples for SIM and NFW with only two 
cases excluded for the NFW model because of a larger virial 
mass. This is a consequence of the fact that the SIM mod- 
els results generally fitted a range of masses (see Table [l]) 
that, according to the Eq.((5]), correspond to 7 ^ 1, i.e. more 
"cuspy" systems. 

On the contrary, almost all the lenses enter the G sam- 
ple in the Burkert case because they have a far lower virial 
mass and the three missing lenses are excluded since they 
have log Mvir ^11. When one strengthens the constraint on 
the precision of the recovered (o"ap, Me), the SIM model be- 
comes clearly preferred, its B sample being the richest one. 

The SIM mod el turns out to be favored by the AIG 
and BIG statistics (|Liddle II2004I ) which measure the signifi- 
cance of the best-fit taking into account the different num- 
ber of degrees of freedom by penalizing the introduction of 
unnecessary parameters. These two estimators are defined 
respectively as : 

yl/C= -21nA^ax + 2n , 

BIC = -21n£max +nlnAfdata , 

where £max is the maximum likelihood value, n the num- 
ber of model parameters and A'data the number of con- 
straints. Assuming a Gaussian likelihood (as we have implic- 
itly done), it is — 21n£max = Xmim while n = 1(2) for the 
SIM (NFW and Burkert) model with A^data = 2. The SIM 
model turns out to have the lowest AIC and BIC values and 
then the favoured with respect the NFW and Burkert ones, 
although the small values of AAIC = AICsim — AICnfw 
and ABIC = BICsim - BICnfw makes the NFW model 
statistically equivalent in most cases (with the Burkert 
model generally excluded because of the highest values). 

A possible caveat is in order here. Although formally 
defined in the same way, the values for the SIM, NFW 
and Burkert models, are not fully equivalent since, for each 

^ Actually, tlie concentration definition strictly applies to tiie 
NFW case. We decided to re-define it for the Burkert model for 
convenience, although we do not assume that this has a defined 
physical meaning. 



lens, the total errors (eapjes) depend also on the adopted 
model because of the propagation of the uncertainties on the 
stellar masses and effective radius (see GTIO for detailf[f|). 
Although such differences are actually quite small and do not 
impact significantly the relative ranking of the models, we 
have nevertheless computed the percentage best fit residuals 
and found (for the good samples lenses) : 

rms [{afp' - a^) /a°l'] = 7% , 13% , 14% ; 
rms [{M^"' - M'^) /A/^""] = 12% , 19% , 16% ; 

for SIM, NFW and Burkert models, respectively. Again, we 
find that the SIM model best reproduces the observed data 
so that we can safely argue that it should be preferred with 
respect to the NFW and Burkert models. This is one of the 
central results of this paper, since we seem to have proved 
that a mono-parametric density profile (including most of 
the physics of the galaxy collapse) works better than two- 
parameter density profiles with a much weaker physical con- 
tent (e.g. collapse and evolution of coUisionless DM parti- 
cles). However, this result is based on a different procedure 
of fitting, thus one can argue that the final significance of 
the fit might have been affected by the use of binned data to 
constraint the halo parameters. In the next Section we will 
use the well known correlations between the two density pa- 
rameters of NFW (and Burkert) profiles to check whether 
their one-parameter re-writing can match the one-by-one 
data with similar significance. 

3.1.2 SIM vs one-parameter NFW and Burkert profiles 

For the NFW density profile there is a well established cor- 
relation between the virial mass and the concentration as 
found in N - body simulations. Following the recent analysis 
in Mufioz- Cuartas et al. (2011), we therefore set : 

log cnfw = aNFw(2:) log (^J~T^J[^ + fcNFw(z) (9) 
with 

aNFw(2) = 0.0292 - 0.097 , 

, , 110.001 2469.720 

'^^"^("^ ^ -2-1- 16.885 + (2 + 16.885)2 ' 

for 10 ^ log [Mvir/(ft"^ Mq)] < 15. If we neglect the scatter 
in Eq.@ as a first approximation, this one parameter ver- 
sion of the NFW model can be fitted to the SLAGS lenses 
with the virial mass as only unknown quantity. Using this 
approach, we thus find a viable solution for all the 59 ob- 
jects in the sample with 53 (50) lenses satisfying the selection 
criteria used to define the good (best) sample for the SIM 
case. In particular, the quality of the fit may be quantified 
by noting that 

rms(l - ~ 4% rms{l - M^V^D - 13% 

^ Summarizing, the difi'erence between the theoretically predicted 
and observed values of the velocity dispersion and projected mass 
entering the evaluation are normalized with respect to the er- 
rors, but these errors are partially model dependent. For instance, 
the uncertainty on the theoretically predicted Mp^oj is propa- 
gated from the one on the stellar mass through a formula which 
depends on the adopted DM halo profile. 
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Figure 1. Best fit virial mass (left) and inner slope (right) for tiie SIM model with Salpeter {x axis) and Chabrier IMF (y axis). 



for the good sample and 

rms(l - a%/a:l') ~ 4% rms{l - M^^ /M^l ~ 11% 

for the best one. We have repeated the same exercise 
with the Burkert model, by using relation between the core 
density po and the core radius ro found in Salucci & Burk- 
ert (2000) to reduce the density model to one-parameter 
profile. Imposing this relation and fitting each single lens 
using the virial mass as only parameter, we have found a vi- 
able solution for only 6 lenses, which suggests that a mono- 
parametric Burkert profile is uncapable to match the SLAGS 
data and we do not discuss this case hereafteiQ. 

If we go trough a more detailed comparison of the SIM 
to the NFW model, we have seen that the latter works better 
in matching the aperture velocity dispersion, while the SIM 
is in a slightly better agreement with the projected mass 
values. Overall, however, both models well match the data 
so that the choice of which model is most viable should be 
driven by physical motivations. We indeed prefer the SIM 
model since the mass— central slope relation comes out from 
a physical model rather than being the outcome of numerical 
simulation. However, this argument alone does not allow us 
to definitely abandon other density models and as far as 
both the SIM and one parameter NFW (as well as the two- 
parameter version of NFW and Burkert models, as discussed 
above) well reproduce the SLAGS lens data, we will hereafter 
consider all options and a posteriori check the difference 
introduced by the different models. 

3.2 Changing from Salpeter to Chabrier IMF 

As fiducial values for the stellar mass of each lens, we have 
adopted the estimates given in Table 4 of Auger et al. (2010) 
under the assumption of a Salpeter IMF. However, since the 
IMF is still an uncertain variable in the mass analysis, in this 
Section we want to investigating the impact of a different 
IMF choice on the SIM results. 

We have therefore consider the case of the Ghabrier 
(2001) IMF which returns stellar masses smaller by a factor 
1.8 with respect the ones obtained with the Salpeter IMF. 



" Note, however, that this is likely a consequence of having as- 
sumed that the po - ''0 relation found at 2 = applies to any 
z. 



As the Ghabrier IMF provides the lowest masses compatible 
with the colors of the galaxy this will allow us to minimize 
the contribution to the velocity dispersion and projected 
mass by the stellar component. 

As for the Salpeter IMF, it turns out that all the lenses 
may be well fitted with the best fit theoretical values of tJap 
and Me within 2a (la) of the observed ones for 59 (53) out 
of 59 lenses in the SLAGS sample. Seemingly, the ability 
of the SIM model to fit the data is not affected by the IMF 
choice. However, when looking at the estimated virial masses 
for the Ghabrier IMF case, the inferred values turned out 
to be generally larger (i.e., logMvir > 14.0) than the ones 
obtained by the Salpeter IMF such that only 25 (24) lenses 
enter the G (B) sample. 

This is made clearer in Fig.[TJ where we show that 
the best fit virial masses for the SIM + Ghabrier model are 
almost one order of magnitude larger than those for the 
SIM -|- Salpeter case because of the lower stellar masses pre- 
dicted by the Ghabrier IMF which maximize the DM contri- 
bution to the central parts. This is "seen" by the SIM model 
as the presence of a more cuspy profile which implies larger 
7 and larger Mvir according to Eq.([5ll. 

Based on the argument of unrealistic virial masses we 
are inclined to rule-out the SIM -I- Ghabrier model. This 
is however becoming a common conclusion from different 
analyses: e.g., gravitational lensing (see, e.g., Treu et al. 
2010, Gardone & Tortora 2010) and studies of the central 
DM fraction in local ETGs (Napolitano et al. 2011) suggest 
that observations are consistent more with a Salpeter rather 
than a Ghabrier IMF for ETGs (unless some strong AG is 
considered). Moreover, it has been also suggested that the 
IMF might vary w ith luminosity (|Renzini fc Giotti I [1993I : 
iTortora et al. l2009l ) with a Salpeter one being generally pre- 
ferred for brighter systems (as probed by the SLAGS lenses) . 



4 DARK MATTER SCALING RELATIONS 

Despite the large uncertainties, we have shown that the SIM 
halo model is able to fit the combination of galaxy kine- 
matics and lensing data fairly well. We have also fixed the 
stellar mass contribution by assuming the Salpeter IMF, we 
can now start looking into the DM properties of the galaxy 
sample. 

The scaling relations we are interested in are commonly 
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Figure 2. Best fit relations between tlie effective column density SdmC^o) and the total luminosity L^, stellar mass M*, effective radius 
Re and halo mass M200 from left to right. Data in sample B are plotted as points, while the other data not present in B, but in samples 
G and F are plotted as triangles and boxes, respectively. We do not plot error bars to not clutter the plot. However, typical errors 
are (o-(log M,)/ log M^) ~ 1%, {cr(log Mvir)/ log Mvir> ^ 2%, ((r(log<SDM(i?o))/ log5DM(iJo)> ^ 11%, while logL^^ and logi?e have a 
negligible error. The continue black line is the best fit made using sample G (i.e. points and triangles). The red and orange dashed lines 
are the best fit obtained in CTIO, using a NFW + Salpeter IMF and a Burkert + Sal peter IMF, resp e ctivel y. The continue blue line and 
cyan shaded region arc the best fit and the region enclosing the data for the results in lTortora et al.H 1I2OIOI ) using a spherical isothermal 
sphere, SIS (adopting a Salpeter IMF and projected quantities), while the dashed ones arc for the same data, but using masses in A+09. 
Continue and dashed gray lines are the best fits of local ETGs in T+09, using a SIS and a constant M/L profile, respectively. 



written as power laws and can be conveniently converted in 
linear relations in a log -log space, logy = ^ogA + Slog a; 
with comparable uncertainties on (log a;, log y). As a best 
fitting procedure, we will follow the approach as in CTIO and 
adopt the Bayesian method described in D'Agostini (2005 
see also Hogg et al. 2010). 

For the scaling relations we want to examine, we need 
to choose a reference radius where the mass quantities are 
evaluated. As these estimates will imply some model extrap- 
olation, this choice is critical sine results can be strongly 
model-dependent . 

Since our data mainly probe the region close to _Rc, 
this seems the natural choice as the reference radius. Some- 
times, also the core radius Rc, introduced in the cored den- 
sity models, and the same Rs, where the logarithmic slope 
of the density profile is —2, are taken as reference radius. It 
is, however, easy to check that both these quantities are far 
larger than _Rc so that one would stay in the inconvenient po- 
sition of deriving column density estimates in regions much 
more far away of the distance constrained by the data ex- 
tension. As the proposed model has proven itself to fairly 
well match the observations around Ra we will compute the 
scaling relations at this reference distance. 

As a final remark, we stress that using Mvir as a model 
parameter allows us to avoid extrapolation of the mass esti- 
mates constrained at Rc out to the far larger virial radius, 
as commonly done in parametric modeling procedures. This 
makes our virial quantities intrinsically more robust. 



4.1 Column density 

We start by considering the correlations of the DM column 
density 5dm (i?) ~ M^^^{R)/-kR^ with luminosity, stellar 
mass, efi'ective radius and virial mass. The best fit relations 
obtained for the SDmiRe) El are: 



° Hereafter, we will discuss only the results for the sample G hav- 
ing checked that fully consistent constraints are obtained using 



log (7?c) = 0.06 log (Lv/10" Lq) -f 3.17 , 

log5DM(i?c3) = 0.08 log (M,/10" Mq) + 3.12 , 

log5DM(i?o) = 0.05 log i?e +3.13 , 

log {R,) = 0.10 log (M200/IO'' Mq) + 3.02 , 

which are also shown in Fig.[2]as compared to the individual 
galaxy datapoints. The intrinsic scatter of the fitting proce- 
dure, (Tint, is formally negligible for all the scaling relations 
above and will not to be considered further on, while the 
true uncertainty of the fit is dominated by the data point 
scatter, an issue we will come back later on. As a virial 
mass estimate we use M200 (i.e., the mass within the radius 
R200 where the mean density is 200 times the cosmological 
mean matter density) for homogenity with previous litera- 
ture studies. 

Taking at face values, these relations suggest a non uni- 
versality of 5dm (J?o) in agreement with what we have found 
in CTIO using a different halo model but the same reference 
radius. A similar comparison with other results in literature 
is not straightforward because of difference in the model 
adopted, the radius where Som{R) is evaluated for each in- 
dividual system, and finally the galaxy sample considered. 
E.g., using a NFW halo profile, Boyarsky et al. (2009, here- 
after BOO) found a strong correlation between 5dm (-Rs) and 
the halo mass over on wide range of masses, ranging from 
dwarf spheroidal galaxies to galaxy clusters. On the con- 
trary, Donato et al. (2009, hereafter D09) using the Burkert 
(1995, B95 hereafter) profile found a remarkably constant 
5dm (J?c) {Rc being the core radius of the B95 model) with 
the luminosity on a sample of local spirals and ellipticals. 
These "cored" profile are in contrast, though, with our SIM 
model which has no core but an inner cusp. Furthermore, as 
shown in CTIO, an inconvenience of the B95 model is that it 
does not allow a decent fit to the SLACS data unless one as- 
sumes unreasonably low values of the virial M/L ratio. Our 

the F and B samples. This can also be qualitatively seen in the 
figures where all the points are plotted with different symbols. 
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68% CL 


95% CL 




(x) 




68% CL 


95% CL 


<5dm(-Rc) - Lv 


0.06 


0.02 


0.02 


(-0.09, 0.14) 


(-0.16, 0.23) 


3.1670 


3.1671 


3.1682 


(3.1631, 3.1694) 


(3.1597, 3.1729) 


5dm(Ro)-M* 


0.08 


0.06 


0.07 


(-0.04, 0.17) 


(-0.16, 0.28) 


3.1204 


3.1228 


3.1157 


(3.0740, 3.1782) 


(2.9801, 3.2478) 


•SdmC^g) - -Re 


0.05 


0.07 


0.07 


(-0.04, 0.19) 


(-0.17, 0.30) 


3.1257 


3.0788 


3.0506 


(3.0117, 3.1882) 


(2.9229, 3.2960) 


oriAvi ( jLo ) - Ajonn 

^ D Wl \ / ^"ZUU 


0.10 


0.10 


0.10 


(0.05, 0.16) 


(-0.02, 0.21) 


3.0189 


2.9892 


2.9615 


(2.9294, 3.0745) 


(2.8660, 3.1560) 


iStot - Lv 


-0.17 


-0.19 


-0.19 


(-0.27, -0.11) 


(-0.37, -0.02) 


3.4390 


3.4344 


3.4364 


(3.4253, 3.4394) 


(3.4245, 3.4424) 


5tot - A'h 


-0.16 


-0.16 


-0.16 


(-0.24, -0.07) 


(-0.34, 0.02) 


3.5390 


3.5207 


3.5144 


(3.4777, 3.5747) 


(3.4353, 3.6331) 


Stot - Re 


-0.37 


-0.34 


-0.33 


(-0.45, -0.24) 


(-0.59, -0.13) 


3.7635 


3.7391 


3.7487 


(3.6724, 3.8000) 


(3.5613, 3.9340) 


iS*-^*- — A//onA 


0.01 


0.01 


0.01 


(-0.06, 0.08) 


(-0.12, 0.16) 


3.4469 


3.3381 


3.3760 


(3.1088, 3.5102) 


(3.1013, 3.6076) 


9dm - Lv 


-0.04 


-0.02 


-0.01 


(-0.17, 0.12) 


(-0.30, 0.28) 


-9.628 


-9.625 


-9.626 


(-9.631, -9.615) 


(-9.635, -9.614) 


goM- Mi, 


-0.03 


0.00 


-0.01 


(-0.14, 0.13) 


(-0.29, 0.29) 


-9.607 


-9.621 


-9.611 


(-9.702, -9.550) 


(-9.777, -9.441) 


9dm - fie 


-0.04 


0.03 


0.03 


(-0.12, 0.18) 


(-0.33, 0.35) 


-9.590 


-9.637 


-9.625 


(-9.771, -9.535) 


(-9.921, -9.332) 


9dm - A/20C) 


0.25 


0.24 


0.24 


(0.17, 0.30) 


(0.10, 0.37) 


-9.947 


-9.978 


-9.961 


(-10.130, -9.864) 


(-10.145, -9.770) 


g*-Lv 


-0.42 


-0.38 


-0.38 


(-0.51, -0.26) 


(-0.63, -0.13) 


-9.394 


-9.395 


-9.395 


(-9.397, -9.394) 


(-9.399, -9.391) 


g*-^L 


-0.35 


-0.28 


-0.28 


(-0.41, -0.14) 


(-0.58, 0.08) 


-9.180 


-9.217 


-9.191 


(-9.287, -9.179) 


(-9.422, -9.066) 


9* -fie 


-1.03 


-1.02 


-1.02 


(-1.12, -0.92) 


(-1.24, -0.80) 


-8.475 


-8.455 


-8.428 


(-8.529, -8.426) 


(-8.647, -8.334) 


9* - A/200 


-0.07 


-0.07 


-0.08 


(-0.15, 0.00) 


(-0.25, 0.10) 


-9.308 


-9.305 


-9.306 


(-9.378, -9.231) 


(-9.527, -9.091) 



Table 2. Constraints on the slope and the zeropoint of the correlations involving <SDM(fie), <Stot, 9dm and 3* using the sample G. For 
each parameter, we report the best fit, mean and median values and the 68 and 95% confidence ranges. Results from the fit to the other 
samples are fully consistent so that we will not report them, but make them available on request. 



lens systems, instead, is made of (mainly) ETGs at interme- 
diate z, thus partially overlapping with B09 and completely 
complementary to the D09 sample. For these systems, we 
cannot exclude that some of the discrepancy with the re- 
sults above can be advocated to some evolution with the 
redshift of the slope of the 5dm (-Re) vs {Lv, -M*, Rc, -M200). 
Albeit not detailed investigated elsewhere (but see, e.g., B09 
and NRTIO for some hints), the galaxy morphology can be 
another important responsible of the discrepancies in the 
scaling correlations above. Our data are therefore not the 
best to solve the differences between B09 and D09 contrast- 
ing results, however they might provide a benchmark result 
for the correlations expected for ETGs under a generalized 
halo model with no fixed cuspy density profile (somehow 
covering the intermediate range between the B95 and the 
NFW central behavior). 

This is particularly true due to the uncertainties on 
the slope and zero-point of the fitted relations. Although 
we use a robust Bayesian fitting method, the confidence re- 
gions around the parameter fit are large (see Table [2]) due 
to the large error bars on the DM related quantities. Such 
uncertainties make any assessment on the non-universality 
of the relations rather weak. Indeed, considering the 68% 
confidence ranges, a zero value for the slope can be statisti- 
cally excluded only for the Sdm{Rc)-M2oo relation, which 
is however consistent with zero within 95% CL. 

Thus, looking at the correlations with the stellar quanti- 
ties (-Re, Mt, Lv), the SoMiRe) universality cannot be ruled 
out. On the other hand, the strong correlation with the halo 
mass -M200 in Fig.[2]works against the Som{Rc) universality. 
This might be one reason for the contrasting results between 
D09 and B09: in fact the former show the absence of a trend 
of the column density with the luminosity, while the latter 
show a strong trend with the halo mass, which is along the 
same line of our conclusions. 



However, even if we look at the correlations with lumi- 
nosity, we are inclined to argue that the claim of the con- 
stant Sdm{Rc) with luminosity is motivated for the later 
type galaxies, while the ETGs have an intrinsically larger 
scatter (see also NRTIO) which is the effect of a stronger 
correlation with the halo mass (see also tj4.2p . 

Focusing on the difference produced by the adoption 
of the different halo profiles, we can now compare the re- 
sults obtained with the SIM model with findings in recent 
literature using the same IMF and SLAGS dataset. 

GTIO have found a lower average DM column densities 
over the sample, A log 5dm (i?e) ~ 0.1 — 0.3 with respect the 
one obtained with the SIM. This is due to the fact that the 
SIM model has an average central slope, 7 ~ 1.0 — 1.3 which 
is steeper than the one of the NFW and implies a projected 
DM mass within Rc larger than the one obtained with the 
same NFW. 

However, the slope of the Sdm{Rc) vs {Lv , Mi,, M200) 
correlations stay very similar for the two halo models (see 
solid black and red dashed lines in Fig.[2]lf| which means that 
the details of the central DM density are possibly insensible 
to the global quantities. On the contrary, the Sdm{Rc) seems 
much more sensitive to the scale of the luminous matter as 
the slope of the 5DM(^?e)-^Ze correlation is inverted with 
respect the NFW. In particular, from Fig. 2, we see that 
the NFW produces almost similar 5dm (-Re) for smaller -Re 

^ We note that in particular the correlation with M200 might 
look somehow too shallow to an eyeball check. We have checked 
that if we use the sample F, the best fit slope increases to 0.13 
which is closer to the 0.17 value found when the NFW model is 
used (red dashed) similarly to the correlation found with Lv and 
Mi,. Moreover, a steeper slope could be found by excluding the 
least massive points, although there seems to be nothing unusual 
for these lenses to be reasonably excluded here and not elsewhere. 
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(red dashed line) which means that the two density profiles 
are very similar for more compact objects, while it clearly 
produces much smaller DM column density at larger _Re, 
where the SIM model accommodates the steeper slopes. De- 
spite the low significance of the trend, this is a key result 
of our analysis because this points to a correlation between 
the halo central slope and the size of the luminous matter in 
the sense that bigger galaxies are formed in more cuspy ha- 
los. This is q ualitatively consistent with w hat already found 
m DP09 and lDel Popolo fc Kroupal (|2009l ). which have sug- 
gested that more massive haloes are cuspier. As a side note, 
the presence of steeper cusps for massive systems goes in 
the same direction of the effect expected for the canonical 
AC which seems to be necessary to model extended kine- 
matics, e.g., us ing planetary nebulae as a mass tracers, in 
massive ETGs ("Napolitano ct al |201ll ) and not in more reg- 
ular ones ((Nap ohtano ct al. ,20o3)- If so, the steeper cusps 
would reduce the strength of the actual AC recipe with 
respect the classica l prescriptions jBlumenthal et al~lll986l : 
ICnedin et al. I [20041 ) as reported elsewhere. 

Going to a comparison with more general mass profiles, 
T-l-10 adopted a singular isothermal model to describe the 
stellar + DM mass profile. Here the difference is that the 
stellar masses have been derived using different population 
models. The solid blue lines in the two middle panels of 
Fig-El shows that their SumiRa) - t^'h and St>m{Rc)-Rc re- 
lations have the opposite trend of the one we have found 
from the SIM. Note that a decreasing iSDM(-Re) with both 
the stellar mass and the effective radius is also found by 
T-l-09 for local ETGs, adopting the same mass model and 
stellar population properties (gray continue line in Fig. [2]), 
but fitting very central velocity dispersions. A similar (al- 
beit shallower) decreasing trend and a lower average column 
densities (of A log >SDM(i?e) ~ 0.3 — 0.5) are also found from 
T-f 10 for a model with a constant M/L model (gray dashed 
line in Fig. [2]). 

As T-f 10 do not assume a halo density profile, but de- 
rive the DM mass by subtraction of the stellar component 
from the total density profile, their approach seem much 
more dependent on the stellar population analysis (i.e., IMF, 
priors on stellar population parameters, availability of pho- 
tometric or spectral data, etc.). E.g., while stellar masses in 
T-l-10 are, on average, consistent with the ones from A-l-09, 
some tilt is present and finally propagates to a different final 
estimate of iSdm(-Rc). Indeed, adopting the same total mass 
model, but the A-l-09 masses, the SumiRe) - A'h correlations 
changes its sign, while the SnmiRe) - Re one becomes shal- 
lower so that the disagreement with our SIM based results 
is partially alleviateq^. This check is a warning to us about 
the importance of the mass model and stellar population 
analysis adopted which, in principle, is an important player 
in the definition of the DM scaling relations. 

Keeping this in mind, we decided to check the scaling 
relations of the total column density 5tot (obtained including 
the stellar mass within Re) obtained with the SIM model for 



^'^ Note that the results obtained in T-flO, using the two stellar 
mass sets, populate similar regions in the spaee cSDM(^e)-Af* 
and SuM{Re)-Rc (see cyan region and the one enclosed within 
the blue lines in Fig. [2]). 



our particular stellar population and IMF choice. As best fit 
correlations we find : 

log5tot = -0.17 log (Lv/10" Lq) -f 3.44 , 

log5tot = -0.16 log (A/*/10" Mo) + 3.54 , 

log5tot = -0.37 log i?e + 3.74 , 

log5tot = 0.01 log (Maoo/lO'^ Mq) +3.34 , 

whose marginalized constraints on the slopes and the zero- 
points are given in Tabled A zero slope is excluded at the 
68% CL for <Stot vs (Lv,M*,iic) correlations, while we do 
not find any correlation with M2oo- Thus, the net effect of 
the inclusion of the stellar component is to tilt the trend al- 
most uniformly clockwise toward more negative slopes (and 
contemporary increase the zeropoint due to the addition of 
the stellar mass). In fact, we have checked that the slope in 
the central regions of our light profiles become shallower at 
larger _Re, possibly combining with the DM slope the results 
above. One of the possible driver of such trends is galaxy 
merging, which has been shown to produce shallower slopes 
in simulated haloes (e.g.. iBovlan-Kolchin fc Ma|[20o3 '). The 
trend of the Stot s eems consistent with the recent finding 
from lAuger et all (|2010l ') which, adopting a simple power- 
law density to shape the full mass profile of the SLAGS 
lenses, have found a strong inverse correlation of the total 
density slope with Re- 

4.2 Newtonian acceleration 

The issue of the universality of the column density S-dm{R) 
has some important dynamical consequences if one considers 
that it can be easily related to the Newtonian acceleration 
gDuir) = GMT3M{r)/r^ (Gentile et al. 2009, hereafter G09). 
In particular, G09 have shown evidences for the universality 
of both g-DM(Rc) and gt,{Rc) with the label DM {*) referring 
to DM (stellar) quantities, over a (small) sample of spirals 
and ellipticals. Let us first discuss the DM case evaluating 
the acceleration at the effective radius to be consistent with 
our choice throughout the paper. For the best fit relations 
(plotted in Fig.O , we get : 

loggDM(-Re) = -0.04 log (Lv/lO" Lq) - 9.63 , 

loggomiRe) = -0.03 log (Af*/10" Mq) - 9.61 , 

loggDM(i?e) = -0.04 log iio - 9.60 , 

\ogguM{Re) = 0.25 log (M200/IO'' Mq) - 9.95 , 

where the accelerations are in units of mjs^. Similarly 
to the column densities, we find that the DM Newtonian 
acceleration has no correlation with the stellar quantities 
{Lv , Mi, Re), while it strongly correlates with M2oo- Thus 
we confirm that goM is a constant with respect to the stellar 
quantities although it is larger than the value found by G09 
for the goM computed at R^. The latter is an upper limit 
to our Re estimates since generally Rc ^Rc and for a given 
constant core density po the Newtonian acceleration scales 
linearly with radius, i.e. guM{Rc) ^ 5dm(^?c). As already 
noted by NRTIO, this result works against the universality 
of the goM, which instead seems to scale with the morpho- 
logical type (and possibly the mass). 
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Figure 3. Best fit relations between the DM Newtonian acceleration (in units of m/s^) and the total luminosity Ly, stellar mass M*, 
effective radius Rc and halo mass Af200 from left to right. Typical error bars are (cflog gDM(Ke)]/ log gDM(J2e)) — 2%. Symbols and lines 
are the same as in Fig. [2] 



The SIM estimates obtained here are discrepant with 
the resuhs from CTIO for the NFW and B95 (Fig.[3ll. In 
particular our goM are located in between the two reference 
halo models, while the difference in the slopes with all the 
quantities are statistically insignificant. This result is some- 
how different with respect the SoyiiRc) in Fig. 2 where the 
CTIO estimates are almost everywhere lower than the SIM 
estimates, which might be tracked to the fact that here we 
are considered 3D quantities. In particular, the large val- 
ues obtained by CTIO for the B95 are related to the fact 
that typical core radius in CTIO are smaller than Rc thus 
not representing a real cored profile, but rather more re- 
sembling a pseudo-isothermal sphere, with a rather steep 
slope around Rc- In fact the deviations are more marked 
for the smaller systems which shall have also shallower in- 
ner slopes (according to Eq. 5), and thus are overestimated 
by the steeper densities implied by the B95 models. On the 
other hand, the gDM(J?o) are systematically larger than the 
values obtained in CTIO for NFW because of the steeper 
3D slopes. We remark here that the 3D quantities are more 
weakly depending on the inner slopes with respect the pro- 
jected ones shown in Fig. 1, a fact that translates with an 
null correlation of the ^dm with _Rc. 

Going to the stellar Newtonian acceleration, we find as 
best fit relations : 

log5,(i?e) = -0.42 log (Lv/10" Lq) -9.39 , 

log5*(i?c) = -0.35 log (M*/10" Mq) -9.18 , 

log5*(i?o) = -1.03 log i?o -8.48 , 

log5„(i?o) = -0.07 log (M200/IO'' Mo) -9.31 , 

which clearly demonstrate that this is not a universal quan- 
tity (and indeed the confidence ranges for the slope in Table 
[2] exclude a zero slope at the 68% CL) in agreement with 
CTIO. Note, however, that a correlation of (;*(i?c) with the 
stellar quantities (Ly,M*,_Ro) is expected, being all these 
quantities involved in its definition. On the contrary, the cor- 
relation with M200 is very weak, which shows that there is 
not any strong dependence of the stellar mass in the centers 
on the global DM content. 



4.3 Dark matter mass content 

In the previous sections we have seen that the total and DM 
column densities correlate differently with luminosity, stellar 
mass and _Rc which shall depend on the different spatial 
distributions of stellar and dark matter. Here we want to 
investigate the consequences of the column density scaling 
relations in term of the central DM fractions. To this end, 
we note that 

'^DM ^ M'^^{Rc) ^ fP"'°-if_R ) 

Stot M.P^°'iR,) + M^Z'iR^) 

where /dm^(^o) is the projected DM fraction. Considering 
the best fit values of the correlations of both column den- 
sities with stellar parameters, we obtain that the projected 
DM fraction within R^ scales with both the luminosity, stel- 
lar mass and i?o, as f^Z^i^") oc L^^^, f^Z^iR<=) « M-,°-^'^ 
and/P'j^j(7?e)oci?e"-^2^ 

The same results are found by direct fitting of the same 
quantities of the individual galaxies, adopting our Bayesian 
fitting procedure : 

log/P'^°J(J?e) = 0.21 log (Lv/lO" Le) - 0.21 , 
log/SM^(J?c) = 0.18 log (A/./lo" Mo) - 0.32 , 
log/gM^(i?c) = 0.41 log i?e - 0.58 , 

while the 68% confidence ranges of the slopes are 
(0.15,0.26), (0.12,0.24) and (0.34,0.45) for the f^Z'-^v, 
/dm"'"-^* and fjjM-Re, respectively. We have therefore a 
clear evidence that the different scalings of Som{Rc) and 
5tot with the stellar quantities are an expected consequence 
of the varying DM content within the effective radius. 

Similarly, we define the three dimensional DM fraction 
as /dm(-Rc) = Mdm(-Rc)/[M*(_Rc) + Mom{Rc)] and investi- 
gate its correlation with the stellar quantities. Our best fit 
relations are : 

log/DM(i?o) = 0.22 log (Lv/lO" Lo) - 0.44 , 
log /DM(i?o) =0.15 log (M*/10" Mo) -0.53 , 
log/DM(i?c) = 0.59 log J?e - 0.98 , 

while the 68% CL for the slope are (0.08, 0.30), (0.02, 0.24), 
(0.36, 0.71) respectively. Despite the large uncertainties, 
these results show that brighter, more massive and bigger 
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systems have a larger DM content within the effective ra- 
dius in qualitative agreement with previous results in liter- 
ature, regardle ss the adoption of deproje c ted or projected 
quan tities (e.g. Padmanabhan et al71l2004l : [Cappellari et alj 
I2OO6I : T-H09: lAuger et al.ll2010l : NRTIO: T+10). We, how- 
ever, note that our best fits relation with Ly and M* are 
much shallower than what is found in CTIO for the fiducial 
NFW model ( and Salpeter IMF), bu t in good agreement 
with the ones in lCardone et alH (|2009l ). where a phenomeno- 
logically motivated general halo profile was used to fit the 
same SLACS data considered here. 



4.4 Impact of the IMF choice 

All the results discussed so far has been obtained using 
the SIM model for the dark halo and a Salpeter IMF. 
As discussed in Sect. 3.3, choosing a Chabrier IMF leads 
to SIM models having a virial mass larger than 10^* Mq 
for more than 50% of the sample so that we have pre- 
ferred to exclude the SIM + Chabrier combination when dis- 
cussing the scaling relations. This is a reasonable choice 
under the hypothesis of a universal IMF which does not 
depend on galaxy parameters. There are different evi- 
dences that this might not be the case, although there 
is still not any consensus whether the IMF might change 
with galaxy morpholog y, luminosity / mass and/or stellar 
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population paramet e rs (iDave I 20081: van Dokkum 
Holden et aTl I2OI0I: iNapolit ano. Romanowskv fc Tortora I 



2010l : lvan Dokkum fc Conrovll2010l . I2OIII '). Here we decided 
to check what can be the impact of the IMF choice, by 
computing the scaling relations for a mixed sample made 
out of the 25 lenses in the G sample of the SIM -I- Chabrier 
model and the remaining 26 lenses of the G sample of the 
SIM -|- Salpeter model. In a sense, we are here postulating a 
IMF varying with the halo virial mass and approximating 
such a variation with a rough step function. 

Let us consider first the DM column density within R^. 
For the best fit relations, we get : 

log5DM(i?e) = -0.003 log (Lv/lO" Lq)H-3.21 , 

log5DM(i?e) = -0.02 log (A/*/10" Mq) +3.22 , 

log 5dm (-Re) = -0.09 log i?e + 2.33 , 

log {Re) = 0.06 log (A/200/10'' Mq) + 3.11 . 

Comparing with the values in Table 1, we see that the 
best fit relations are shallower for (L^jAf*) and steeper 
for (-Re,Af20o). However, if we consider the 68% confidence 
ranges, the slopes are fully consistent so that the change in 
the slope can not be considered statistically significant. 

Similarly, the total column density iStot turned 
out to have the best fit relations with slopes 
(-0.21,0.05,-0.04,-0.04) for {Ly, M-,, R^, M200) re- 
spectively. Also for these relations the confidence ranges 
significantly overlap with respect to the Salpeter IMF (see 
Table 1) thus the two cases do not differ significantly. 

For the DM Newtonian acceleration we obtain : 

log5DM(i?e) = -0.05 log (Lv/10" Lq) -9.45 , 
log5DM(i?e) = -0.07 log (M,/10" Mq) -9.42 , 
loggDM(i?e) = -0.13 log i?e -9.33 , 



log 5DM(-Re) =0.13 log (M200/IO'' Mq)-9.67 

which are comparable with the results for slope and zero- 
point as in Table 1. This suggests that the universality (or 
lack of) of the DM Newtonian acceleration is not signifi- 
cantly affected by the choice of the IMF (a similar results 
is found for the stellar Newtonian acceleration which we do 
not report for brevity). 

Finally, we have considered the DM mass fraction 
within Re and found: 

log /DM(iie) =0.18 log (Ly/lO" Lq) -0.30 , 
log/DM(i?e) = -0.10 log (A/,/10" Mq) - 0.24 , 
log/DM(i?e) = 0.44 log i?e - 0.70 . 

As expected, these results turned out to be significantly 
different from the Salpeter IMF case. In particular, the /dm - 
Lv and /dm - Re relations turned out to be shallower, while 
the /dm - has a negative slope and a significant anti- 
correlation. This is the consequence of the step IMF func- 
tion assumption that produced a significant increase of the 
/dm (-Re) for the lower mass systems with the Chabrier with 
respect the smaller fuM{Re) obtained with the Salpeter IMF 
assumed for the more massive ones. The same argument ap- 
plies to the /dm - Lv and /dm - Re relations which turned 
out to be shallower. 

As a final remark, we can conclude that scaling rela- 
tions are overall mildly affected by the IMF assumption in 
the dynamical/lensing analysis, with correlations being com- 
parable among results obtained either assuming a universal 
(Salpeter) IMF, or a non universal mass dependent IMF 
(Chabrier for less massive systems and Salpeter for more 
massive ones). A universal Chabrier IMF seems to be ruled 
out because it produces too many systems with unrealisti- 
cally large virial masses for galaxy systems. 



5 CONCLUSIONS 

While there is a general consensus on the ubiquitous pres- 
ence of DM in galaxies, there is still an open debate about 
its mass density distribution is and its role it has exactly 
played in the galaxy formation scenario. These are crucial 
issues since the exact density distribution and the assembly 
processes that dark halos have undergone through might tell 
more on the actual nature of the DM itself. Scaling relations 
among DM related quantities and stellar properties may give 
important hints about the relative interplay between the two 
main constituents of galaxies. 

In this context, it seems particular intriguing the exis- 
tence of some universal properties of the DM quantities like 
the column density or even a presence of a characteristic 
accel eration scale (e.g. iDonato et al. 1 l2009l : iGentile et al. I 
which might be related to common formation pro- 
cesses on many mass scales or morphological categories. Re- 
cent works have argued either in favor or against the pres- 
ence of such universal values (B09, CTIO, NRTIO). Dealing 
with DM properties, though, means necessarily to deal with 
indirect, model dependent quantities, which can make the 
conclusions on these parameters strongly affected by the 
particular model adopted. One approach might be to use 
models which bracket the widest range of DM properties de- 
rived from the N-body cosmological simulations, e.g. from 
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the cuspy profiles predicted by tlie classical NFW to the 
"cored" profiles of the B95 (see e.g. CTIO). The disadvan- 
tage of this approach is that it gives a partial vision of how 
the DM scaling relations vary in the two extreme regime 
without allowing to generalize the results in case the actual 
DM properties are in the between. 

In the attempt of testing more general, theoretical mo- 
tivated DM halo profiles, like the one proposed by the sec- 
ondary infall model as implemented in DP09, we have mod- 
eled the central velocity dispersion and the projected mass 
within the Einstein radius of a large sample of ETGs lenses 
at intermediate redshift {{z) ~ 0.2). 

In particular, the DP09 models adopted a modified sec- 
ondary infall scenario including (in a semianalytic way) the 
effect of angular momentum, dynamical friction and adia- 
batic collapse of baryons. As a first important result, we 
have shown that this model is fully compatible with obser- 
vations of ETG and well performing in the data fitting. This 
is a significant step forward with respect to previous analy- 
ses where SLAGS sample have been modeled with standard 
NFW (CTIO), as the SIM model has the main property to be 
fully assigned by a single parameter (the virial mass) instead 
of the two parameters required by more "standard" halo 
density models like NFW (which is specified once concen- 
tration and virial mass are given) or Burkert profile (char- 
acterized by the core radius and the central density). 

Even if we consider the corr elations between the 
halo parameters for the NFW |Navarro et al.~l 1 19961 : 
iBullock et"al~ll200ll : [Hennawi et al. IbOOTl ) and Burkert pro- 
files (e.g. Salucci & Burkert 2000) which allow us to re-write 
these models as a function of one parameter only (e.g. the 
virial mass), we have seen that the best-fit to the observed 
quantities turned out to be generally poorer than the one 
provided by the SIM model (see Sect. [3]), with only NFW 
providing somehow similar significance than the SIM model. 
This is mainly because in the typical mass range spanned 
by the ETG sample considered here, the SIM model predicts 
cuspier profiles according with Eq.((S}, i.e. NFW-like or even 
cuspier. 

With this novel DM setup, we have estimated the 
DM column density SomiRe), the Newtonian acceleration 
SDM(-Re) and, finally the DM mass fraction fuM{Re) along 
with their correlations with the stellar total luminosity, mass 
and size. The best fit relations show that SoM{Ra) is almost 
constant (possibly increasing with, if any) over the range 
of stellar parameters {Lv,M^,,Rc) probed by the adopted 
dataset, while it is clearly strongly correlated with the halo 
mass M200. This result is actually consistent with both G09, 
claiming a characteristic density scale over a large range of 
galajcy luminosities, and with B09 which instead has found 
a correlation with the dark halo mass. 

Similarly we have found an even more remark- 
ably constant DM Newtonian acceleration gDM{Rc) with 
{Lv , Mi,, Re), and still a strong correlation of this quan- 
tity with the virial mass. In this case, though, the abso- 
lute constancy of the gr>M{Rc) does not allow to justify the 
correlation with the virial mass but probably says more of 
some intrinsic properties of the dark matter halos, and of 
the non-existence of an universal acceleration scale. In fact, 
we have confirmed here a former evidence from NRTIO, that 
ffDM(i?e) of the ETGs is on average larger than the one ob- 
tained for late-type galaxies (as e.g. in Donato et al. 2009, 



having considered that their column density obtained at the 
core radius are an upper limit for the same quantity if com- 
puted at the Re). 

As such, one could argue of the existence of a Newto- 
nian acceleration growing with the morphological type and 
(possibly) with the stellar mass (if including all the Hubble 
sequence) as well as shown by the trend with Al200- 

We need to conclude this reasoning with a caveat, un- 
derlying all analyses based on model dependent approaches: 
all DM quantities correlations against (Ly, Mi,, Re) are crit- 
ically dependent on the adopted halo model and stellar IMF 
(in our case a Salpeter IMF) . A direct comparison with pre- 
vious literature results is complicated by systematic effects 
due to differences in both the halo model (SIM vs NFW 
or Burkert) and the radius where the DM quantities are 
evaluated (-Rc vs Rs or Rc). In particular, the comparison 
of results obtained with the SIM model with similar works 
based on the adoption of the more standard NFW and B95 
has highlighted an interesting correlation between the cus- 
pyness of the dark halo and the size of the parent galaxies. 
In the SIM approach, in fact, larger galaxies seem to assem- 
ble in dark halo having steeper cusps, while more compact 
massive ETGs are accommodated on shallower cusps. This 
is an interesting hint that might be cross-checked on hydro- 
dynamical simulations and seems to be a crucial test for the 
hierarchical model as a whole. 

Finally we stress that there are still some weaknesses in 
the analysis proposed, mainly posed by the limited galaxy 
sample (e.g. the SLAGS sample only probes a quite limited 
range in both stellar luminosity and mass), and the large 
uncertainties on the derived quantities. However we have 
considered this a benchmark test for more general dark halo 
model laws based on a simple physically motivated galaxy 
model. This is the first step of a broader plan to improve the 
testing along different roads. On one hand, a detailed the- 
oretical investigation is needed to find out quantities which 
depends as less as possible on both the adopted halo pro- 
file and stellar IMF. Similarly, one should check whether the 
choice of Re as a reference radius where DM quantities are 
evaluated is the most convenient one finding a compromise 
between the need to not extrapolate outside regions directly 
probed by data and the halo model characteristics. On the 
other hand, stronger constraints on the slope of the investi- 
gated correlations could be obtained by narrowing the un- 
certainties on the DM quantities. The use of a one parameter 
model, like the SIM, is expected to help by eliminating the 
degeneracies among halo parameters which generally plague 
the analysis adopting the NFW profiles or similar, and con- 
tribute to the overall error budget. A further improvement 
would be obtained fitting the full velocity dispersion profile 
rather than its aperture value only, although such a strategy 
could be applied to local ETGs only. Finally, a larger sample 
spanning a wider range in [Ly , Mi,, Re, M200) would allow 
to further narrow down the confidence ranges for the slopes 
of the investigated correlations by both improving statistics 
and better tracking the different trends. 

Should both constancy with stellar luminosity, mass 
and size and independence on the fitting procedure details 
be successfully demonstrated, one could safely conclude that 
a proposed quantity is indeed universal and use such a result 
to constrain galaxy formation and evolution scenarios. 
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